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This thesis presents some results of the nonlinear dynamics experiments in an argon dc glow 
discharge plasma. The thesis has been divided into six Sections. 

Section 1 contains an introduction to nonlinear dynamics experiments in plasmas, and the time 
series analysis. The study of nonlinear dynamics or chaos theory has emerged in the last three 
decades or so as an important interdisciplinary area of research encompassing a wide range of fields 
like: fluids, plasmas, biomedical sciences, finance, turbulence, astronomy, material sciences, etc. 
The advent of fast computing systems have been of great help in solving nonlinear equations which 
were intractable analytically. Low dimensional deterministic chaos has become a part and parcel of 
almost every field of engineering like, mechanical, electrical, civil, chemical, etc., and this has lead 
to the development of new tools and data analysis techniques, like nonlinear time series analysis, 
wavelet transform, etc. 

In an electron beam plasma system in 1984, Boswell showed that the natural oscillations of an 
electron beam propagating parallel to a magnetic field, went through a period doubling bifurcation 
to a chaotic state with increase in the beam current. In 1989, the intermittency route to a chaos 
was observed in the low frequency self oscillations in the undriven DC discharge system. Homoclinic 
chaos was reported in the same electrical discharge plasma system in which the deterministic chaos 
had been observed by T. Braun, et al. Period adding route to chaos and period subtracting had been 
observed in on ion beam plasma in a double plasma device. Other nonlinear phenomena like mode 
locking, period pulling, etc., had been observed by Klinger et al. In many other experiments almost 
a similar phenomena had been observed where different types of gases, geometric configurations and 
parametric regimes were explored. 

In all the above experiments, the common feature was that the plasma dynamics goes towards a 
chaotic state as the discharge current or voltage (DV) was raised. On the other hand, the inverse 
route, i.e., chaos to ordered state transition has not been observed experimentally in plasmas. The 
possibility of such a route in a plasma sheath had been shown by Tomejiro Yamagishi and Makoto 
Tsukayama [Tomejiro Yamagishi and Makoto Tsukayama, J. Phys. Soc. Jpn. 69, 2883 (2000)]. 
numerically, in the presence of dissipation. Generally, with dissipation or decrease in the input 
energy, the system goes to an ordered state. In our experiments the chaotic behavior of the floating 
potential fluctuations, goes to a regular state with increasing DV. 

In Section 2, we have described the experimental setup, diagnostics, and data analysis procedures. 
The experiments were preformed in an argon dc glow discharge plasma. Plasma was produced inside 
a cylindrical hollow cathode dc glow discharge system of radius 45 mm and cylindrical anode of radius 
« 1 mm placed inside the cathode concentrically. The entire electrode system was housed inside a 
vacuum chamber whose pressure (p) could be varied between 0.001 to 1 bar. A dc discharge voltage 
applied between the two electrodes produced the plasma, and a Langmuir probe located between 
the anode and the cathode was used to monitor the floating potential fluctuations, and a black 
and white CCD camera was used to measure the size of the anode glow. The noise (HP33120A) 
and a signal generators (Fluke PM5138A) were coupled with the discharge power supply for the 
noise induced experiments. The measured plasma density and the electron temperature were about 
10 7 cm -3 and 3 — 4 eV respectively, and the corresponding electron plasma frequency (f pe ) and the 
ion plasma frequency (f P i) were about 28 MHz and 105 kHz respectively. The analysis of the floating 
potential fluctuations has been carried out using the linear and nonlinear time series analysis. In the 
linear analysis, power spectrum, autocorrelation, probability density function, etc., and in nonlinear 
analysis Lyapunov exponent, correlation dimension, Hurst exponent, etc., have been estimated. 
Some statistical tools like 'normal variance' and 'absolute mean difference' have also been used to 
characterize coherence and stochastic resonances. Most of the data analysis were carried out using 
the software MATLAB. 

In Section 3, we have presented anode glow related observation of chaos to order transition 
and homoclinic bifurcation in the glow discharge plasma. Depending upon the gas pressure (p), 



a discharge was struck at different discharge voltages (DV), which exhibits a Paschen curve like 
behavior. The region to the left of the Paschen minimum is narrower than the region to the right 
and the behavior of the plasma floating potential fluctuations was different on these two sides of the 
Paschen curve. To right of the Paschen minimum, an anode glow or anode spot was observed to 
form on the anode, whose size decreased with the increase in the DV, before finally disappearing at 
a certain DV. The annular radii of the glow estimated from the experimental observations was in the 
range of f .3 — 0.32 mm. This is consistent with the theoretical estimations of thickness (S) using the 
relation S « 3.7 x 10^ 6 -^p, where k, T, at, and P are the Boltzmann constant, room temperature 

in Kelvin scale, ionization cross-section in cmT 2 and pressure in mbar respectively. The estimated 
thickness of the anode glow using the above relation for P « 0.95 mbar, T w 300 K, is w 0.81 mm 
is within the range of the thickness as obtained from the experiment (1.3 — 0.32 mm). 

The floating potential fluctuations have been analyzed for three typical pressures in the region 
greater than the Paschen minimum. At about 0.89 mbar (pd w 20.02 mbar-mm), the discharge 
was initiated at w 288 V and simultaneously irregular relaxation type of oscillations were observed 
but turned into regular oscillations with increase in the DVs. Around 509 V, these oscillations 
disappeared and this point (DV) is termed as the bifurcation point (Vh). A similar behaviour was 
observed at 0.95 mbar (pd — 21.37 mbar mm) and at 1 mbar (pd = 22.5 mbar mm ). Initially, the 
power spectrum of the oscillations are of broad band nature indicating a chaotic behavior of the 
system. An estimate of the frequency of these instabilities from the ion transit time in the plasma 

r(d) — — d/ J '^ L , where d is the electrode distance, is « 19 kHz which agrees quite well with 
the frequency of the relaxation oscillations of the floating potential. 

The relaxation oscillations have been attributed to the formation of highly nonlinear structures 
like double layers, and hence, we tried to estimate the correlation dimension (D corr ) and the +vc 
Lyapunov exponent (\l) of all the signals. The estimated D corr for the three pressures (0.89, 0.95 
and 1.0 mbar) were greater than 3.8 to begin with and decreased with increase in DV. Since D corr is 
a measure of the complexity of the system, initially, the system was in a complex state as D corr for 
all the three pressures were high and decreased with increase in DV.The presence of a +ve Lyapunov 
exponent (Xl) is the most reliable signature of the chaotic dynamics. Positive \l was obtained for 
DV 283, 284, and 290 V at 0.95 mbar and for 293, 296, 300 and 305 V at 1 mbar respectively. At 
higher pressures we find that in general \ L is +ve and D corr > 3, suggesting a low dimensional 
chaos. Wavelet analysis also showed a presence of chaos at the initial stage of the discharge. 

Generally, the floating potential fluctuations were complex in nature at the initial stage of dis- 
charge and became regular with increase in DV and converted to relaxation oscillations before totally 
vanishing at the bifurcation point (Vh). 

The time period (T) of these relaxation oscillations increases with increasing DV, before eventually 
vanishing beyond a critical DV (Vh), and for larger values, the autonomous dynamics exhibits a 
steady state fixed point behavior. It is observed that the autonomous dynamics undergoes an 
exponential slowing down. It is observed that ln|V — Vh| vs Time is a straight line, indicating 
that the system goes through a homoclinic bifurcation. The floating potential fluctuations exhibit 
relaxation oscillations on the one side of the Vh (bifurcation point) and a stable fixed point on the 
other side, called the excitable state and is useful to study noise invoked dynamics— stochastic and 
coherence resonances. 

In Section 4, the noise invoked dynamics has been investigated. When the DV is modulated 
with a stochastic signal, the autonomous dynamics will be recovered and this phenomena is termed 
as coherence resonance. For stochastic resonance, the DV was perturbed by a combination of a 
subthreshold periodic pulse train and a Gaussian white noise. The time series of the system response 
in the presence of an identical subthreshold signal for three different amplitudes of imposed noise. 
It is observed that there is little correspondence between the subthreshold signal and the system 
response for a low noise amplitude. However, there is an excellent correspondence at an intermediate 
noise amplitude. Finally, at higher amplitudes of noise, the subthreshold signal is lost amidst 
stochastic fluctuations of the system response. Absolute mean difference (AMD) has been used to 
quantify the information transfer between the subthreshold signal and the system response and the 
unimodal structure is the signature of the SR phenomena. 

In Sections 5, we have presented the observation self organized criticality behavior in plasma. 
When the system was operated in the pressure region less than the Paschen minimum, for a small 
range of p (0.9 - 1.5 x 10~ 2 Torr or 1.2 - 2 x 10~ 2 mbar), it was observed that the behavior of 
the floating potential fluctuations was consistent with self organized criticality (SOC). In order to 
establish the SOC behavior, we had checked the power law behavior of the power spectrum, and 
the presence of the long-range correlation by estimating Hurst exponent (H) (using R/S technique) 
and the exponent (a) of autocorrelation function (ACF) decay, and the nongaussian probability 
distribution function. The results of the Hurst exponent (0.96 ± 0.01) being greater than 0.5, the 
ACF exponent, a ~ 0.30, nongaussian PDF, and power spectral index (3 « 1.7 ± 0.1 in the pressure 
range 9 x 10~ 3 — 1.6 x 10~ 2 Torr, are consistent with the systems exhibiting SOC like behavior. 
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Finally, we have presented a summary and future plan of work in Section 6. Though we have 
investigated some very important and interesting nonlinear processes during the the course of this 
thesis work, there are, many problems which we plan to take up as a part of our future activities, 
are as follows. 

• Effect of the noise on the autonomous dynamics. 

• Effect of suprathreshold signal and noise to the plasma. 

• Investigation of the existence of canard orbits in plasma. 

• Non-chaotic attractors may be investigated by applying two non-commensurate periodic signals. 

• Chaos control and synchronization. 

• Modeling of the experimental results. 

Dusty plasma which can be produced very easily in our experimental system, is another area 
where, a lot of nonlinear phenomena can be explored. 
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1. OVERVIEW OF NONLINEAR DYNAMICS 
WITH RELEVANCE TO PLASMA PHYSICS 

1.1. Introduction 

Nonlinear dynamics and chaos theory started with the 
intention of investigating the qualitative behavior of non- 
linear problems which were difficult to solve analytically. 
In the first half of the 1900s, scientists were interested 
in nonlinear oscillators like the van der Pol oscillator, 
and with the discovery of the high speed computers 
many different nonlinear problems were solved pQ. In 
1963, Lorenz (5J introduced, the fluid convection model 
to study the dynamics of the atmospheric weather, which 
exhibited chaotic behavior in the numerical simulations 
and this model later came to be known as the famous 
Lorenz model. It shows very complicated dynamics like 
chaos, inverse period doubling bifurcation, etc., depend- 
ing on the model parameters. Chaos theory flourished 
when Reuelle and Takens [20] proposed a new theory 
for the onset of turbulence in the fluids and Gollub and 
Swinney verified it experimentally in a concentric ro- 
tating cylinder using water, and it was observed that 
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with the increase in the rotation rate the radial veloc- 
ity fluctuations went to a chaotic state through a period 
doubling bifurcation, which were analyzed using Fourier 
techniques [21j. In the late 70s, Feigenbaum discovered 
the universal constant called the Feigenbaum constant to 
characterize the universal features of the period doubling 
bifurcation [22J. Later, bifurcation diagrams, Lyapunov 
exponent, correlation dimension, etc., derived on the ba- 
sis of chaos theory, have been used to characterize chaos 
and its different routes (23] El] E5J. Since the develop- 
ment of nonlinear dynamics, most applications have been 
in the field of fluid dynamics, particularly fluid turbu- 
lence and since fluid and plasma turbulence arc closely 
related, the concepts of nonlinear dynamics had been suc- 
cessfully utilized in plasmas. 

In the late sixties and later, the idea of van der Pol 
oscillator was applied to explain the growth and satura- 
tion of the plasma instabilities which could not be done 
using conventional linear theories Ell E3 EE] EH ES3 ISU 
[SlISllMIMJIMllSllMllMllinilllJIllllSlIll]- Abrams 
et ai, first observed the nonlinear phenomena like pe- 
riod pulling, frequency cntrainment, etc., in periodically 
forced self oscillatory plasmas [IS]. B.E. Keen, et ai, 
showed using a two fluid model that the ion sound in- 
stabilities behaved in a manner similar to a van der Pol 
oscillator when subjected to a driving force E5]. By the 
mid eighties researchers began looking for experimental 
evidences of deterministic chaos in plasmas. Experimen- 
tally, the period doubling bifurcation and chaotic behav- 
ior were observed in an electron beam plasma system 
in 1984, wherein, Boswell showed that the natural os- 
cillations on an electron beam propagating parallel to a 
magnetic field, went through a period doubling bifurca- 
tion to a chaotic state with increase in the beam cur- 
rent [3J. The period doubling route to chaos was also 
reported in the driven pulsed filament discharge plasma 
and an electrical discharge tube [461 Hi]- I* 1 1989, the 
intermittency route to a chaos was observed in the low 
frequency self oscillations in the undriven DC discharge 
system [4] [5] . Homoclinic chaos [6] was reported in the 
same electrical discharge plasma system in which the de- 
terministic chaos had been observed by T. Braun, et 
al, Quasiperiodic route to chaos was observed in the 
driven and undriven discharge plasmas 5, 48 1 49] in mag- 
netic multipole devices [15] US]- Period adding route to 
chaos and period subtracting had been observed in ion 
beam plasma in double plasma device [7]. Other non- 
linear phenomena like mode locking, period pulling, fre- 
quency entrainments etc., had been observed by Klinger 
et al H |3_5] ESI S3 [30] ED 133 [33]. In arc plasma chaos 
have been observed by S. Ghorui, et al., [M] [33] 135] , 
In many other experiments almost similar phenomena 
had been observed where different types of gases, geo- 
metric configurations and parametric regimes were ex- 
plored [9] HO] El H3 H31- I n all the above experiments 
the common feature was that the plasma dynamics goes 
towards a chaotic state as the discharge current or dis- 
charge voltage (DV) was raised. 



On the other hand, inverse route, i.e., chaos to ordered 
state transition has not been observed experimentally in 
plasmas. Possibility of such a route in plasma sheath 
had been shown by Tomejiro Yamagishi and Makoto 
Tsukayama [14] numerically, where dissipation was im- 
portant to lead to such situations. Generally, with dissi- 
pation or decrease in the input energy, the system goes to 
an ordered state. In our experiments the chaotic behav- 
ior of the floating potential fluctuations generated during 
the creation of an anode glow, goes to a regular state with 
increasing the discharge voltage. This happens because 
with increase in DV, the glow vanishes and hence the sys- 
tem becomes stable and this has been discussed in detail 
in Section [3] 

With the advent of nonlinear dynamics, many other 
interesting phenomena like controlling chaos [57 , noise 
induced resonances 58J, synchronization of chaotic sys- 
tem |59j . quantum chaos, etc., have been observed. Some 
of these phenomena are also observed in plasma, specially 
coherence resonance near Hopf bifurcation was studied by 
Lin I and Jcng-Mci Liu [60 and A. Dinklagc, et al. |61j. 
Our experiments have been very conducive for investiga- 
tion both stochastic and coherence resonance phenomena 
and have been discussed in Section [4] We are presenting 
results of both types of resonances in the glow discharge 
plasmas near homoclinic bifurcation in Section [4] 



1.2. Time series analysis 

The instabilities that we have studied in this thesis 
were recorded in the form of time dependent floating po- 
tential fluctuations. To understand the characteristics of 
the instabilities, one has to analyze the fluctuations using 
the time series analysis which has been given here. 

A time series is a sequence of data points of an ob- 
served variable at equally spaced time intervals and time 
series analysis comprises of methods that attempt to 
understand such time series, often either to understand 
the underlying context of the data points (where did they 
come from? what generated them?), or to make fore- 
casts (predictions) . Applications cover virtually all areas 
of statistics, but some of the most important ones in- 
clude economic and financial time series, analyzing and 
modeling and forecasting of physical time series (hydrol- 
ogy, earth sciences, astronomy, oceanography and biol- 
ogy, astrophysics, plasma, etc.), simulation of time series 
models and examination of their properties, filtering and 
smoothing of any type of time series, statistical estima- 
tion of parameters of the time series models, either linear 
or nonlinear, with several types of optimization methods, 
etc. 

The objectives of a time series analysis are diverse, 
depending on the background of applications. Statisti- 
cians whose fundamental task is to unveil the probability 
law that governs the observed time series, usually view 
them as a realization from a stochastic process, With 
such a probability law, we can understand the under- 
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lying dynamics, forecast future events, and control them 
via suitable intervention. Time series analysis can be cat- 
egorized in two classes: linear and nonlinear time series 
analysis. Methods for linear time series analysis are of- 
ten divided into two classes: frequency-domain and time- 
domain methods. The former are spectral and wavelet 
analysis while the latter consists of estimating the auto- 
correlation and cross-correlation in the time series. 

Linear analysis tools treat all the irregular or chaotic 
behavior as a stochastic processes though these instabil- 
ities may be generated by deterministic dynamics. Very 
often, we also encounter such systems which have no 
model and and one has to study such a system from a 
single observable. So to overcome this lacking feature 
of the linear analysis and to analyze complex system, 
one needs resort nonlinear dynamics approach especially 
chaos theory. The most direct link between chaos theory 
and the real world is the analysis of the time series from 
real systems in terms of nonlinear dynamics [62] . Based 
on chaos theory, nonlinear tools have been developed to 
characterize such complicated data. The scope of these 
methods ranges from invariants such as Lyapunov ex- 
ponents, correlation dimensions, etc. [231 HH H5] > which 
yield an accurate description of the structure of the sys- 
tem to statistical tools which allow for classification and 
diagnosis even in situations where determinism is not ob- 
vious. 

We have used both types of analysis techniques (linear 
and nonlinear) to analyze the time series data obtained 
from our experiments, because some features that could 
not be revealed by one type of technique was possible 
by the other. A detailed discussion on the data analysis 
techniques has been presented in Section [23] of Section [2] 



EXPERIMENT, DIAGNOSTICS, AND DATA 
ANALYSIS PROCEDURES 



2.1. Experiment 

The experiments were carried out in a glow discharge 
argon plasma produced by a dc discharge in a cylindrical 
hollow cathode electrode system with a typical density 
and temperature ~ 10 7 and 2 — 6 eV respectively. The 
schematic diagram of the electrode system is shown in 
Fig. [T] and its detail is presented in Subsection |2.1| The 
electrode assembly was housed inside a vacuum cham- 
ber as shown in Fig. [2] The neutral pressure inside the 
vacuum chamber was controlled by a needle valve and 
the range of the gas pressure in these experiments were 
between 0.001 to a few mbar. The discharge voltage 
(0 — 999 V) was applied between the cathode and the an- 
ode keeping the anode grounded. A noise (HP33120A) 
and a signal generator (Fluke PM5138A) were coupled 
with the power supply for the noise induced experiments. 
The main diagnostics was the Langmuir probe (described 
in Subsection 2.2 1 used to obtain the floating potential 



(DV) and pressure were the control parameters. Fig. [3] 
is the laboratory view of the whole experimental system 
and the details of the various parts are presented below. 
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FIG. 1: Schematic diagram of the cylindrical electrode sys- 
tem of the glow discharge plasma. The probe is placed at a 
distance I ~ 12.5 mm from the anode. 

By varying the neutral pressure (p), the breakdown of 
the gas occurs at different DV as shown in Fig [3] The 
breakdown voltage (Vt, r ) initially decreases with increase 
in pd, (where, d is the radius of the cathode) and then be- 
gins to increase with pd after going through a minimum 
value resembling a typical Paschen curve |16[ 1631 . This 
curve classifies the operating regions of p and DV for the 
experiments presented in this thesis. When the discharge 
has been operated for pressures greater than the Paschen 
minimum, the floating potential fluctuations show chaos 
to an ordered state transition and has been discussed in 
Section ?? and for less than the Paschen minimum, fluc- 
tuations show self organized criticality behavior and has 
been discussed in Section [6] Chaos to an ordered state 
transition is related to the formation of an anode glow 
or an anode spot on the anode [TH1 [53, 65 . Evolution of 
the glow was monitored by a CCD camera. The ordered 
state, i.e., when the floating potential is in a steady state 
which is basically an excitable state is very useful to per- 
form noise invoked resonance experiments and has been 
presented in Section [4] 

The detail descriptions of the different parts of the ex- 
perimental setup have been given in the following sub- 
sections. 



Electrode system 

The experiment has been carried out in a hollow cath- 
ode cylindrical electrode system as shown in Fig [5] and 
schematic diagram in Fig. [I] Fig [5] (A) and [5] (B) are 
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FIG. 2: Schematic diagram of the whole experimental setup. 1. probe connection to oscilloscope, 2. view port, 3. power 
supply, 4. Langmuir probe, 5. anode, 6. cathode, 7. suspension of electrode system, 8. upper flange, 9. needle valve, 10. 
bottom flange, 11. SS304 bellow, connection to the rotary pump, 12. connection to pirani gauge, 13. toughened glass inside 
the view port, 14. vacuum chamber, and 15. argon gas cylinder; NG: noise generator, and SG: signal generator. 



the side and front view of the electrode system respec- 
tively. The outer stainless steel (SS) cylinder [Fig 5^c)) 
is the cathode and the inner SS cylindrical rod placed 
concentrically with the cathode [Fig |5^d)] is the anode. 
Cathode was kept covered with the nylon [Fig[5je)]. Elec- 
trical connections have been given to the electrodes by 
two connectors shown in Figs[5jal) and (a2). Almost all 
the experiments have been carried out using a cathode 
of diameter 45 mm. 



Power supply 



Discharge voltage has been supplied to the electrode 
system by a dc power supply [Figs. [2] and [3]. It is a 
variable power supply of range 0-999 V and maximum 
output current is 1 amp. Over current protection and 
short circuit protection are available. 



Vacuum Chamber 



SS304 vacuum chamber consists of three parts, bottom 
flange, main cylindrical body and upper flange [Fig. [2]. 
The diameter of the upper and bottom flanges are 300 
mm each and height of 300 mm. We have used two upper 
flanges. It consists of two NW16 coupling through which 
a gas inlet valve and a Pirani gauge have been connected 
and few NW25 coupling for inserting electrode system, 
probes etc. The bottom flange has a pumping port. The 
whole chamber is installed on an MS stand [Fig.[3]. There 
are four side ports each of diameter 200 mm. These ports 
may be covered using toughened glass or blank flanges or 
other suitable flanges fitted with diagnostics. 



Pumping System 



The pumping system is a mechanical rotary pump pro- 
vided by Hind High Vacuum, Bangalore. This pump has 
been mounted on the bed of a MS table [Fig. [3] . The vac- 
uum chamber is connected to the rotary pump through 
an SS bellow followed by a butterfly valve. During pump- 
ing the butterfly valve is kept open, and closed when the 
system is not under operation. The Pirani gauge head is 
connected on the top flange of the chamber to measure 
the pressure. 
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FIG. 3: The experimental setup: 1. Langmuir probe, 2. pirani gauge, 3. gas inlet valve, 4. view port, 5. signal generator, 6. 
noise generator, 7. signal and noise amplifiers, 8. amplifier power supply, 9. argon gas cylinder, 10. discharge power supply, 
11. Langmuir probe supply, 12. pressure meter, and 13. rotary pump. 
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FIG. 4: Vbr vs pd (Paschen curve) for our experimental sys- 
tem. The minimum occurs at (1.69 mbar-mm, 251 V). 
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FIG. 5: Experimental electrode system: A. side view and B. 
front view of the surface; al and a2- connection to the anode 
and cathode, c. cathode, d. anode, and e. nylon coat on the 
electrode. 



2.2. Diagnostics 

Langmuir probe 

A Langmuir probe is a conductor introduced in the 
plasma to measure density, floating potential fluctua- 
tions, etc. It consists of a cylindrical tungsten wire of 
diameter 0.5 mm of length ss 20 mm [Fig. |]. A teflon 
coated stainless steel wire has been soldered to the tung- 
sten wire, and it has been fitted inside a glass tube as 
shown in Fig. [6j The exposed tungsten wire of 2 mm 



length, outside the glass tube is the effective measuring 
area of the probe. It has been inserted in the plasma 
through a side port by NW25 coupling and has the facil- 
ity of axial movement. When it is inserted in the plasma, 
the electrons and ions fall on it. The imbalance between 
the ion and electron current, due to the fact that the 
electrons have a higher mobility between the two, re- 
sults in an accumulation of the electrons on the surface 
of the probe, setting up an electric field. This electric 
field repels the electrons and attracts ions towards the 
probe so that the net current is zero and hence the sur- 
face acquires a negative potential, which is called floating 
potential (</>/). In the present experiments, the floating 
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potential and its fluctuations have been measured using 
a tektronix digital oscilloscope. 




0123456789 10 
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FIG. 6: Typical glass covered Langmuir probe. A. Tip of 
the probe (tungsten wire), B. teflon coated wire, C. soldering 
point of the probe with teflon coated wire, and D. glass cover. 
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Image sensor: 


1/3 B/W CCC 


Picture element: 


EIA: 768X494, CCIR: 752X582 (HXV) 


Sync. system: 


Internal 


o/lN nation. 


iviore tnan 4o ars 


Horizontal Resolution: 


duu J. v iinea 


Min Illumination: 


min 0,1 Lux (F 1.2) 


Auto Gain control: 


Built in. 


Lens mount: 


c/cs 


Video output: 


1 V p.p./75ohm 


Frame Rate: 


30 FPS-NTFC 



TABLE I: Specifications of the CCD camera used in the 
experiments. 



CCD Camera 

Ionization and excitation are both important in glow 
discharge plasmas. The excitation, followed by de- 
excitations with the emission of radiation, is responsi- 
ble for the appearance of glow. Depending upon differ- 
ent plasma parameters, the size of the glow changes and 
hence the evolution of it, is an useful observable for the 
study of the nonlinear behavior of the glow discharge 
plasma. The evolution of the glow has been monitored 
using a black & white CCD camera shown in Fig [7] Size 
of the anode glow has been estimated from the pixel size 
of the digital image taken by this camera. Detailed spec- 
ifications of the camera have been given in Table [TJ 
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FIG. 7: Black and white CCD camera used in the experiment, 
specifications have been given in Table [I] 



2.3. Data analysis techniques 

The analysis of the floating potential fluctuations has 
been carried out using the linear and nonlinear time series 
analysis. In linear analysis, power spectrum, autocorre- 
lation, probability density function, etc., and in nonlin- 
ear analysis Lyapunov exponent, correlation dimension, 



Hurst exponent, etc., have been used. Some statistical 
tools like 'normal variance' and 'absolute mean differ- 
ence' have also been used to characterize coherence and 
stochastic resonances. The analysis has been carried out 
in time as well as frequency domain. Some basics of time 
series analysis techniques and the procedures followed for 
computation have been presented in the following subsec- 
tions. 



Power spectral method 

A time series which has a certain periodicity may be 
represented as a superposition of periodic components 
of sines and cosines. The determination of their relative 
strengths is basically known as spectral or Fourier 
analysis. When the time series is periodic, the spectrum 
may be expressed as a linear combination of oscillations 
whose frequencies are integer multiples of the basic fre- 
quency which is the Fourier series. In order to find the 
characteristic modes, and the dominant mode present 
in the floating potential fluctuations, power spectral 
method has been used. It can also be used to explore 
other characteristic features present in a time series, 
such as, self organized criticality behavior, wherein, 
power spectra shows a power law behavior and it has 
been discussed in chapter ??. Power spectrum has been 
calculated using fast Fourier transform method using 
MATLAB, which has been given below. 



close all 
clear all 
clc 

load data (v) 

(power_v,frequency_v) = power(v,sampling -frequency) 
% power_v and frequency _v are the power and frequency 
% of the signal v which has taken at a sampling rate 
s ampling .frequency 
plot(frequency_v,power_v ) 

function [pyy,f]=power_sig(v,smp_freq) 

% 
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% [pyy, f]=power_sig(v, smpjreq) calculates power 
spectrum 

% of the signal V, sampled at frequency 'smp_freq'. 

%'pyy' is the output 

% power and 'f ' is the frequency. 

% Written Md. Nurujjaman 2006. 

% 
sig=v; 
z=smp Jreq; 
lnth_sig=length(sig) ; 
q=floor(log(lnth_sig)./log(2)); 
y = fft(sig,2<0; 
Py = y* conj(y) /(2«); 
f =z*(0:(2')/2)/(2«); 
pyy=py(l:(2«)/2+l); 



Autocorrelation 

The autocorrelation time of a fluctuating signal is a 
measure of the temporal coherence and is obtained from 
the autocorrelation of the signal. For a time series of 
length n, X = [Xi, i — 1,2, ...n], the autocorrelation 
function (ACF) can be written as 



C(r) 



j + T 



X){X 3 -X) 



(i) 



where X, and r are the mean, and time lag of the time 
series respectively. If there is long-range time dependence 
in the signal, then the algebraic decay of the ACF can be 
written as [57] 



C(r) 



(2) 



for large t, where < a < 1. 

Presence of power law decay of the tail of the ACF 
indicates the presence of long range correlation in the 
time series that is one of the characteristic behavior of 
the self organized criticality [Section [6] . Autocorrelation 
also gives the good estimate of the time delay of embed- 
ding in the phase space that will be discussed in the next 
sections. It has been calculated using MATLAB in-built 
function 'xcorr'. 



Phase space plot 

In order to observe nonlinear structures present in a 
time series and to estimate the correlation dimension, 

the 



Lyapunov exponent, etc. [subsections 2.3 and 2.3 



time series has to be described in a suitable phase space. 
The phase space description provides a powerful tool for 
describing the behavior of a time series in a geometric 
form. As in the real experiments, the observable is the 



sequence of scalar measurements, and one need to recon- 
struct phase space vectors from these scalar observable. 
This can be done by method of delay. 

From a time series (a^, x%, Xjy), where N is the total 
number of points, the m dimensional vector in the phase 
space can be constructed by delay embedding [2S] EH1 EH 
1701171] 



[Xi, Xi+i, a;(i+( m -l)r)J 



(3) 



The time difference (r) between two consecutive com- 
ponents of the delay vectors [Eqn. [3] is referred to as the 
lag or time delay. The minimum embedding dimension 
(m) should be greater than 2Dp for good reconstruction 
of delay vector, where, Dp is the box counting dimension 
of the attractor. A good estimate of r is more difficult 
to obtain. As in the mathematical sense, embeddings 
with the same m but with different r are equivalent for 
noise free data, but in the real case the reconstructed 
dynamics depends upon r. There are several rules to 
choose a time delay. A good time delay may be chosen as 
t = 1 — 1/e, of the autocorrelation of the initial value [62 . 
Then the reconstructed trajectory of the actual dynam- 
ics can be written as X = (Xi; X 2 ; X 3 ; Xm), where 
M = N — (to — 1)t. This trajectory is useful to estimate 
the correlation dimension and Lyapunov exponents that 
have been discussed in the next subsections. 



Correlation dimension 

The correlation dimension of a time series is defined as 
the dimensionality of the space occupied by the points of 
that time series. It can be estimated from the correlation 
integral of the reconstructed trajectory in phase space 
[subsection ??] of the time series. The correlation integral 
can be computed using following equation 



N 



C(r) 



jtEE e(r-MQ-^-l) (4) 



N(N-l) ^ ^ 

v I i=l j=i+l 



where, r is scale length, and is the Heaviside step func- 
tion. Xi, Xj will be obtained using Equation [3] 

The scaling of the function C (r, to) can be written as 

C(r,m)=r D (5) 

where D is the Correlation dimension defined as 

9C(r, to) 



D = lim lim 



OAf^oo d In ? 



(6) 



and can be obtained from the slope of In C(r) vs In r plot. 
We have used the algorithm written by Rosenstein in C 
for computing correlation dimensions. 
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Lyapunov exponent 

Chaotic dynamical systems are sensitive to initial con- 
ditions, and exhibit an exponential divergence in the 
phase space. The divergence can be quantified by an 
exponent called Lyapunov exponent. There are several 
methods to estimate Lyapunov exponent. Here Rosen- 
stein method has been used, because this method is very 
useful for short data length. 

If we consider two points on two nearby trajectories 
of a chaotic attractor, in the phase space, assuming the 
distance between them to be d(0), and after time t, if the 
distance between the two trajectories becomes d(t), then 
the divergence (separation after time t) can be written 
as PS] 

d(t) = d(0)e XLt (7) 

where Xl is the largest Lyapunov exponent. Since a prac- 
tical time series is basically a scalar measurement, the 
largest Xl can be estimated from the reconstructed tra- 
jectories in the phase space. Now if Xj and are the j th 
pair of the nearest neighbor on two nearby reconstructed 
trajectories of that time series in the phase space and dis- 
tance between them is dj(0) = minXj\\Xj — X*\\, then 
the separation after time iAt, where At is the sampling 
time, can be written as 

dj{t) =dj(0)e x ^ iA V (8) 

where dj(0) is the initial separation. The above equation 
can be written as 

In dj (*) = In dj (0) + X L (iAt) (9) 

and the largest Xl can be calculated from the slope of 
the average line defined by 

< \ndj(t) >=< lndj(0) > +X L (iAt) (10) 

Xl quantifies the chaotic nature of the fluctuation, and 
+Ve Xl indicates the presence of chaos. It has been 
calculated using the algorithm developed by Roscnstein 
inC. 



Surrogate method 

The interest in detecting determinism or stochasticity 
in time series arose out of the need to discriminate a 
chaotic time series from a random one. To do this one 
has to first generate surrogate time series from the origi- 
nal time series using some specific rules which is termed 
as surrogate method. Once the property under investiga- 
tion (for instance nonlinearity, determinism, and so on) 
is defined, a method is sought to generate the surrogate 
data that preserves a specific characteristic of the orig- 
inal time series. A discriminant factor is then used to 



compare the original time series with the time series gen- 
erated by a suitably chosen surrogate method. If there 
are significant difference between the discriminating fac- 
tors obtained from the two time series, then one can say 
that the original is not random. 

Several algorithms have been proposed to generate 
the surrogate data in the literature. The most com- 
monly used methods are: Random Shuffled, phase shuf- 
fled and Amplitude Adjusted Fourier shuffled meth- 
ods [73] EH [75] . We have used phase shuffled method 
for our experimental data. The correlation dimension has 
been used as the discriminating factor. The phases of the 
original time series have been randomized by Phase Shuf- 
fled method, by shuffling the fourier phases [73 ] 175 ] 176], 
This method preserves the power spectrum (linear struc- 
ture), but the nonlinear structures are destroyed [75] . 
Now if the original time series has some kinds of non- 
linearity, the correlation dimension estimated from the 
both times series will differ significantly and hence useful 
to detect determinism present in the time series. The 
algorithm written by Michael Small using MATLAB has 
been used to generate the surrogate data and the codes 
may be found in his web page. 



Hurst exponent 

The rescaled-range statistics (R/S) method was pro- 
posed by Hurst and well established by Mandelbrot, and 
Wallis [77 . For the time series defined above [subsec- 
tion ??], the R/S is defined as [77] the ratio of the maxi- 
mal range of the integrated signal normalized to the stan- 
dard deviation: 

R(n) _ max(0, W ll W 2 , W n ) - mm(0, W u W 2 , ...,W„) 

_ (11) 

Here Wk — Xi+x 2 +x 3 + ...+Xk~ kX(n), where X, S 2 (n), 
and n are respectively the mean, variance, and time lag 
of the signal. The expected value of R/S scales like cn H 
as n — > oo, where H is called the Hurst exponent. For 
random data H=0.5, while H > 0.5 for the data with 
long range correlations. H < 0.5 indicates the presence 
of long-range anti-correlation in the data. 

Presence of long range correlation in a time series can 
be estimated using Hurst exponent. The algorithm for 
estimating Hurst exponent was written in MATLAB. 



close all 
clear all 
clc 

load data(v) 

(ln_n,ln_rsl)=hurst_exp(v) ; 
plot(ln_lag,ln_rs); 

% Slope of the fit will give Hurst exponent 
function [ln_n,ln_rsl]=hurst_exp(vl) 

% 
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% [ln_n,ln_rsl]=hurst_exp(vl) Calculates HURST expo- 
nent 

% of a signal vl. %ln_n is the log of lag and ln_rsl is the 

logof(R/S). 

% vl is the time series. 

% exponent of a vector V. 

% Last Modified 02/11/2007 

% MD. NURUJJAMAN 

%=============================== 

v=vl; 
ll=length(v); 

nl=floor(ll/100); % One can define the block smallest 
data length using this. 

i=i; 

while i<=nl 

vi=[]; 
j=i; 

lLf=floor(ll./i); 

% 

while j<=i 

vl(j)=(max(cumsum(v((llJ*(j-l)+l):(ll_f*j))... 
-nanmean ( v ( (11 J* (j- 1 ) + 1 ) : (lLf *j ) ) ) ) ) . . . 
-min (cumsum (v( (lLf* (j- 1 ) +1) : (ll_f*j ) ) . .. 
-nanmean(v( (11 J* (j-1) +1) : (11 _f*j) ) ))))... 
/nanstd(v((llJ*(j-l)+l):(UJ*j))); 

j=j+i; 

end 

% 

rs(i)=nanmean(vl); 

In_a(nl+l-i)=log(U/i); 

ln_rs 1 (nl+l-i) =log(rs (i) ) ; 

i=i+l; 

end 



Normal variance 

Normal variance (NV) has been used to quantify 
the regularity in the signal recovered by introducing 
noise in the plasma. NV [75] is defined as NV — 
std(tp)/mean(t p ) where, std(t p ) the standard deviation 
of the peak to peak distance (t p ) of the recovered signal. 
If the distribution of the peaks are regular NV is small 
and hence it is the estimation of the coherency present 
in the system. NV has been calculated using algorithm 
written in MATLAB. 

function stat=NV(v, delta, smpleT) 

% 

% NV estimates Normalized variance of time series v. 
% NV is defined as NV=standard deviation(v)/mean 
(v) 

% delta is the height of the maximum from its preceded 
value. 

% smpleT is the sampling time of the signal. 

% 



(maxtab, mintab)=peakdet(v, delta); 

% 

% maxJndex = maxtab(:,l); % index of the peaks 
st=smplcT; 

max _indexl=max_index(2:end) ; 
max _index2=max_index(l:end-l) ; 
Max JDiff =max_indexl-max_index2 ; 
P2Ptm=Max_Diff.*st; 
=stat=std(* 2Ptm) /mean(P2Ptm) ; 
% 

function [maxtab, mintab]=peakdet(v, delta) 

% 

%PEAKDET Detect peaks in a vector 

% [MAXTAB, MINTAB] = PEAKDET(V, DELTA) 

finds the local 

% maxima and minima ( "peaks" ) in the vector V. 

% A point is considered a maximum peak if it has the 

maximal 

% value, and was preceded (to the left) by a value lower 

by 

% DELTA. MAXTAB and MINTAB consists of two 
columns. Column 1 

% contains indices in V, and column 2 the found values. 

% Eli Billauer, 3.4.05 (Explicitly not copyrighted). 

% This function is released to the public domain; Any 

use is allowed. 

maxtab = []; 
mintab = []; 

v = v(:); % Just in case this wasn't a proper vec- 
tor 

if (length(delta(:)))>l 

error('Input argument DELTA must be a scalar'); 
end 

if delta <= 

error('Input argument DELTA must be positive'); 
end 

mn = Inf; mx = -Inf; 

mnpos = NaN; mxpos = NaN; 

lookformax = 1; 

for i=l:length(v) 
this = v(i); 

if this > mx, mx = this; mxpos = i; end 
if this < mn, mn — this; mnpos = i; end 

if lookformax 

if this < mx-delta 

maxtab = [maxtab ; mxpos mx]; 

mn = this; mnpos = i; 

lookformax = 0; 

end 
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else 

if this > mn+delta 

mintab = [mintab ; mnpos mil] 

mx = this; mxpos = i; 

lookformax = 1; 

end 

end 

end 



AMD 



Absolute mean difference (AMD) |T7j is the statistical 
tool, proposed to quantify the stochastic resonance in a 
plasma subjected to noise and a periodic signal. AMD is 
defined as AMD = abs(mean(t p /5—l)), and gives AMD 
gives the degree of mimicking the output to the input 
subthreshold signal. Usually, regularity in the stochastic 
resonance is quantified by calculating cross-correlation 
(Co = | < [(arl— < xl >)(x2— < x2 >)] > |) between 
the output and input signal. But in case of plasma it 
is not suitable, because there is always a lag between 
periodic signal that is applied to the plasma and the out- 
put. This lag also varies with time because the plasma 
conditions keep fluctuating over time. Therefore, cross- 
correlation is not the right quantity to be estimated. So 
we have proposed a statistical tool AMD which will be 
independent of lag and is estimated as follows: 

1. First calculate the mean inter-peak distance (5) of the 
periodic signal. 

2. Calculate the inter-peak distances (t p ) of the output 
signal. 

3. Calculate the ((t p — S)/8) 

4. Take the absolute, i.e., | < (t p /8 - 1) > |. 

Therefore, AMD = \ < (t p /5 - 1) > |. 
The algoritm for AMD is almost identical of Normal 
variance. 



3. CHAOS TO ORDER AND HOMOCLINIC 
BIFURCATION IN GLOW DISCHARGE PLASMA 

3.1. Introduction 

In the first chapter, we had mentioned that in almost 
all the chaos related plasma experiments, plasma fluctua- 
tions were observed to go to a chaotic state with external 
control parameters, like discharge voltage, currents, etc. 
On the other hand, there has been no experimental evi- 
dence of inverse routes, i.e., chaos to order. It has only 
been seen numerically in plasma by Tomejiro Yamagishi 
and Makoto Tsukayama [14 . In the present experiment, 
we have reported chaos to an ordered state transition 
for a particular range of plasma parameters given by the 



Paschen curve which has been discussed in Section [2~T1 of 
Section [2] For those pressure regions, with the initiation 
of discharge, a bright glow or spot [Fig [8] was observed 
to form on the anode that was unstable initially and was 
the source of the chaotic fluctuations . So to understand 
the dynamics of the anode glow one needs to analyze the 
floating potential fluctuations that has been discussed in 
next Subsection |3.2| Nonlinear analysis of the fluctuation 
has been carried by estimating the correlation dimension 
and the Lyapunov exponent. 



3.2. Anode glow 

Anode glow is a bright structure or spot of higher 
plasma density, producing anode double layer by addi- 
tional excitation and ionization process and forms around 
an anode or any point on the anode or on a positively bi- 
ased probe [6H E3 [79l [80l EU [82l [83l EH Ell [86l [87]. It 
appears when the discharge current is too low to sustain 
the discharge [64, 65 . Its size and shape differ depending 
upon the plasma conditions. It is related to highly non- 
linear phenomena involving the physics of sheaths, dou- 
ble layers, ionization, beam and possibly external circuit 
interactions. 
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FIG. 8: Evolution of glow size of the anode glow with increas- 
ingDV (a) - (h). 

Figs 8ja) shows that the glow with largest size, ap- 
pears when the discharge is struck at a typical pressure 



13 



of 0.95 mbar and its size decreases with increase in the 
DV until it finally disappears [Figs pta)-|8fti)]. From 
the CCD image analysis the annular radius of the glow 
around the anode, was estimated to be ~ 1.3 mm at 
the beginning of the discharge [Fig |8ja)] and reduced to 
« 0.32 mm [Fig gf)]. 

Theoretically one can obtain an expression for the 
thickness (S) following Bin Song et al. 64J. The rate 
of ion production in the sheath is approximately, 



eV A 



5ai{eV A )N, 



where, n esi m e and N are the electron density, mass and 
neutral density respectively. Va is the positive voltage 
created inside the anode, 5 the thickness of the glow, Ui 
the ionization cross-section depends upon the energy of 
the electron and at 15.76 eV, Oi is « 2 x 10 _18 cto 2 [55] , 
The rate of ion loss from the glow is approximately, 



eV A 



where, rij s , are ion mass and density respectively. 

Balancing the ion loss with the ion production, and 
invoking the condition that rii s ss n es , we obtain 



-6a l (eV A )N « 1 



or 



S « 3.7 x 10" 



. kT 



(12) 



where k, T, and P are the Boltzmann constant, room 
temperature in Kelvin scale, and pressure in mbar re- 
spectively. The estimated thickness of the anode glow 
using the Eqn. [12) for P = 0.95 mbar, T = 300 K, is 
« 0.81 mm which is within the range of the thickness es- 
timated from the image (1.3 — 0.32 mm) shown in Fig [8] 
Fig |9ja) and [9jb) show the tentative model of floating 
potential profile for the present experimental arrange- 
ment. It shows that the floating potential increases 
sharply in the sheath region near the cathode and then 
decreases slowly up to the anode. As the measurement 
near the anode and the cathode was not possible, the 
extrapolated potential profile is shown by dash lines. In 
the the presence of the anode glow, the potential near 
the anode get modified. Experimentally it is observed 
that the anode glow is separated from rest of the plasma 
by a double layer [Ml E3 ES EH] and the height of the 
potential of the core of this glow is of the order of the ion- 
ization potential [89 , which for argon is about 15.76 eV. 
As the anode is grounded, the modified potential profile 
in this case, over the profile shown in Fig[9](a), will be as 
shown in Fig[9jb). As the width of the anode glow de- 
creases with increase in the DVs [Fig [8] the hump shown 
in Fig [9jh) will also shrink accordingly. This shrinkage 
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FIG. 9: Spacial variation of the plasma floating potential: (a) 
without anode glow; (b) with anode glow. Vertical dotted 
and dashed lines are the cathode positions (C) and anode 
position (A) respectively. Horizontal dash dot line is the zero 

level. in the curve near the anode and the cathode are 

the extrapolation of the potential. 



happens because as the electron energy increases with the 
DVs, the necessity of the additional ionization to main- 
tain the discharge decreases [M]- Based on the obser- 
vation of Valentin Pohoat-a, et al. [HI], we feel that the 
relaxation oscillations are double layer or such coherent 
potential structures, which are constantly forming and 
annihilating. 



3.3. Analysis of the floating potential fluctuations 

Chaos to order 

An interesting feature associated with the anode glow 
was the different types of oscillations in the floating po- 
tential at different pressures. We carried out a detailed 
analysis of the fluctuations for three typical pressures 
which are presented here. At about 0.89 mbar (jpd ~ 
20.02 mbar-mm), the discharge was initiated at 288 
V and an anode glow was observed similar to Fig |8ja) . 
Simultaneously irregular relaxation type of oscillations in 
the floating potential were observed [Fig 10 'a)]. Increas- 
ing the DVs, led to an increase both in the amplitude and 
the time period of the oscillations [Figs 10 (b)-|10[i)]. 
The regularity of the oscillation also increases with the 
DVs. Around 509 V both the anode glow and the fluctua- 
tions disappeared simultaneously and this point is termed 
as bifurcation point or critical point (Vh)- At 0.95 mbar 
(pd w 21.37 mbar-mm) the fluctuations are observed to 
be more random [Fig[ll];a)-(i)] than at 0.89 mbar. With 
increasing DV, the final form of the fluctuations before 
their disappearance was the relaxation type of oscilla- 
tion as shown in Figs [TT^Ii) — |TT~|^i) . Increasing the pres- 
sure causes more randomness in the signals as seen in 
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FIG. 10: Sequential change in the raw signal (normalized) at 
0.89 mbar for different voltages: (a) 288 V; (b) 291 V; (c) 295 
V; (d) 301 V; (e) 304 V; (f) 307 V; (g) 327 V; (h) 385 V; (i) 
466 V. All y-axes range form -1 to 1. 
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Fig 12 'a) — (i) at 1.0 mbar (pd w 22.5 mbar-mm) 
served the relaxation type of oscillations at DV 
just before the glow and the fluctuation disappears. So in 
general, the autonomous dynamics shows transition from 
oscillatory to steady state behavior. 
Linear analysis: 

To understand the nature of the fluctuations we have 



analyzed them using Fourier techniques. Figs 13 [I}, {II} 
and {III} show the power spectrum [TH1 [T5] calculated 
for the signals at three different experimental pressures 
shown in Figs [10] [TTJ and [12] respectively. A broad range 
of frequencies (< ion plasma frequency) is observed at 
the initial stage of the discharge, at pressures 0.89 mbar 
[Figfjlri} (a)], 0.95 mbar [FigQlI} (a)] and 1.0 mbar 
[Fig 13 [III} (a) — (b)]. The figures also show that with 
increasing DV, number of peaks present in the power 
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FIG. 11: Sequential change in the raw signal (normalized) at 
0.95 mbar for different voltages: (a) 283 V; (b) 284 V; (c) 286 
V; (d) 288 V; (e) 289 V; (f) 290 V; (g) 291 V; (h) 292 V; (i) 
293 V. All y-axes range form -1 to 1. 



spectrum, decreases. At higher DV, i.e., at the end of 
the sequence of the power spectra we observed one ma- 
jor peak with its harmonic [Figs [l3|[I}(f) — (i) and III(i)] 
or two peaks [Figs 13[II}(i)] just before the fluctuations 
cease. The power spectra clearly show that the system 
goes from an irregular to a more or less regular state with 
increasing DV. Another interesting feature that has been 
observed for all the three pressures is that the dominant 
frequency, i.e., the frequency which has maximum power 
decreases with increase in the DVs [Figs [l4"|a) — (c)]. On 
the contrary, in many experiments [8 , 92, 93J, the domi- 
nant frequency is observed to increase with the DV. The 
basic difference between the present experiment and the 
experiments performed in Ref [H] was in the elec- 
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FIG. 12: Sequential change in the raw signal (normalized) at 
lmbar for different voltages: (a) 293 V; (b) 296 V; (c) 298 V; 
(d) 299 V; (e) 300 V; (f) 305 V; (g) 308 V; (h) 310 V; (i) 312 
V. All y-axes range form -1 to 1. 



trode configuration and biasing arrangement. In those 
experiments [5J [55] IMj the system was planar type is 
which the anode voltage was increased keeping the cath- 
ode at ground whereas ours is a cylindrical system in 
which we increased the negative voltage on the cathode 
keeping anode at ground. In those systems electrons orig- 
inating at cathode get lost at the anode. In our case the 
electrons can undergo multiple reflections and hence in- 
creases the ionization length. The fluctuations are prob- 
ably due to ion acoustic instability driven by electron 
beam plasma interaction and also from electron and ion 
trapping in the potential wells, associated with the anode 
glow [Fig [9]. 

Ions produced inside the anode glow due to collisions 



of the accelerated electrons across the hump with the 
neutrals, makes the glow unstable [65], which is proba- 
bly responsible for the relaxation or the random oscil- 
lations [H [22] • An estimate of the frequency of these 
instabilities can be obtained from t he ion transit time in 
the plasma [8] r(d) = p^— = d/y^jf, where d is the 
electrode distance. The estimated ion transit frequency 
(^) for our experimental system is ~ 19 kHz which agrees 
well with the frequency of the relaxation oscillations of 



the floating potential shown in Figs 13 [I}(f) — (i), {H}(i) 



and {III}(i). The higher frequency oscillations could be 
due to ion acoustic instabilities since the conditions are 
quite conducive to excite these low frequency instabili- 
ties. In the next subsection we have presented nonlinear 
analysis of the fluctuations. 
Nonlinear analysis: 

The presence of relaxation oscillations have been at- 
tributed to the formation of highly nonlinear structures 
like double layers [HI]. We therefore estimated the cor- 
relation dimension (D corr ) and the +ve Lyapunov ex- 
ponent (Xl) °f a ll the signals. In this experiment, it is 
observed that the nature of the fluctuations of the poten- 
tial did not vary for almost the whole day as long as the 
controlling parameters are kept constant, so the analyzed 
signals can conveniently be taken to be stationary. 

The correlation dimension (D corr ) has been calculated 
using well-known Grassberger-Procaccia techniqu es |23l 
172] , and detail has been discussed in Subsection 2.3 of 
Section [2] A typical plot of In C(r, m) vs In r has been 
shown in Fig 15 for the embedding dimension (m), in the 
range of 2 to 10, from which we have estimated the corre- 
lation dimension (D corr ) for a typical signal (Fig 11 'a)) 
at DV w 283 V, at 0.95 mbar. The scaling region had 
to be chosen carefully, since for too small scale lengths 
the correlation sum is heavily distorted by noise and the 
higher scale lengths are limited by attractor dimensions. 
From the above plot the correlation sum exhibits a power 
law behavior within a certain range of r as shown by the 
vertical dotted line. In C(r, m) vs In r plots are almost 
parallel for higher m (i.e. m = 7 — 10) and the corre- 
sponding best fit has been shown by line. D corr vs 

m is also shown in the inset of Fig |15[ It shows that D corr 
becomes constant at higher m and this constant value of 
D corr is the correlation dimension of that particular sig- 
nal and for all of our data D corr > 3.8 to begin with. 
The D corr at 0.89, 0.95 and 1 mbar for different DVs 
have been shown by open circle(o) in Figs [l6|a)-|l6][c) 
respectively. It is observed that there is a decreasing ten- 
dency of Dcorr except for some intermediate values of DV 
at higher pressures, where it is seen to increase and then 
decrease again. Since D cor r is a measure of the complex- 
ity of the system it is likely that the system complexity 
increases at those intermediate DVs. In our experiment 
initially the system is in a complex state as D CO rr for all 
the three pressures is high and decreases with increase in 
DVs [Fig 16 'a) — (c)]. D CO rr ~ 1 just before the system 
reaches stable state finally, is the indicator of the peri- 
odic state and these periodic nature are also prominent 
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FIG. 13: Power spectrum of floating potential fluctuations at initial discharge voltages at filling pressures: (a) 0.89 mbar, (b) 
0.95 mbar and (c)l mbar. 




TABLE II: Error in the estimation of the correlation dimen- 
sion and the +ve Lyapunov exponent 



data length D c 



500 
1000 
1500 
2000 
2500 



9.47 
4.00 
3.85 
3.78 
3.87 



A, 



0.240 
0.056 
0.058 
0.054 
0.058 
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from the raw data and the power spectrum [Figs 10 
12 and Fig [l3|l}— {III}]. D corr shows that the system 
stabilizes itself with increase in DVs. 

The presence of a +ve Lyapunov exponent (A^) is the 
most reliable signature of the chaotic dynamics and it is 



estimated using the Rosestein algorithm [25| which has 
been discussed in Subsection |2.3| of Section [2] Fig. [T7| 
shows < hi dj(t) > vs iAt for m=2 to 10. A clear scal- 
ing region is seen at a higher m shown by line. The 

positive Xl has been identified at 0.95 and 1 mbar for 
some DVs and they are shown in Fig 16 'b) and 16 'c) 



by +ve sign respectively. Figs 16 'b) and (c) show that 
X L becomes positive for 283, 284, and 290 V at 0.95 
mbar and for 293, 296, 300 and 305 V at 1 mbar respec- 
tively. Though D corr quantifies the complexity present in 
the system, it does not guarantee the presence of chaos, 
which is determined by . At 0.89 mbar initially we 
have a high D corr [Fig 16 'a)], but the Xl is not positive 
in this case, implying that the system is not in chaotic 
state. At higher pressures we find that in general Xl is 



17 



5 
4 
3 
2 

1 

288 300 



1(a) 



6-o- 



325 



350 



375 



- i + a 0.052 lis ' 








+ X =0.024 us" 1 
O ' 






"^O _ 


o — 













— 







4 

3 
' 2 
1 

5 
4 
3 
2 

1 

293 295 



(b) 



284 286 288 




290 292 


- + X L 0.01 1 (.is' 1 

+ a. =0.022 us"' ~" 


"*\ 


=0.030 |is"' 


°\ + 






; o^ / 




\ " 

° V 



(c) 



300 305 310 

DV (V) 



FIG. 16: D corr vs discharge voltages. +Ve \l has been shown 
by +ve sign, at (a) 0.89 mbar; (b) 0.95 mbar; (c) 1 mbar. 




15 20 25 30 

time (|is) 



35 40 



FIG. 17: Average In dj (t) for different embedding dimensions. 
The best region for Al has been shown by dotted line. 
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FIG. 18: Surrogate data using phase randomizing method: 
(a) original data; (b) surrogate data; (c) Correlation dimen- 
sions vs m for original (A) and its surrogate data (©). 



We also carried out a wavelet analysis of the fluctua- 
tions at 0.95 mbar and obtained the presence of chaos. 
The ridge plot method (Chandre, et al. [94 ) clearly shows 
that the presence of many ridges [Fig [l9[a)] indicating 
that at the initial stage of DV (283 V) the system is 
chaotic [51], whereas at 284 V [Fig 19 'b)] the system is 
weakly chaotic as only few ridges are present. With in- 
crease in DV the system becomes almost regular in nature 
[Fig 19 b) — (e)] where ridges are approximately constant 



with time. But at 290 V [Fig 19 'f )] abrupt change in fre- 
quency, i.e., presence of transient features suggest that 
these are the effect of nonlincarity present in the sys- 
tem [5S] . From the multifractal analysis [55] we obtained 
the D corr which is shown in Fig. 20 (— • — ). This agrees 
quite well with the D corr [open circle in Fig. 
Grassberger-Procaccia techniques [531 E2] ■ 
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using 



Homoclinic bifurcation 



+ve D corr > 3, suggesting a low dimensional chaos. 

D corr and Al have been estimated for different data 
lengths to check the error on the estimations of the D corr 
and Al- The estimated D corr and Al at different data 
lengths have been shown in Table [ll] Both and A l and 
Dcorr tend to show stable results for higher data lengths 
(larger than 1000) as seen from Table [II] 

In order to validate the above observations, we car- 
ried out the surrogate data analysis using phase shuffled 
method that has been described in Subsection 12 . 3 1 of Sec- 
tion [2j seen in Fig 18 'c) The D corr for the original data 
saturates at higher m, whereas for the surrogate data it 
increases with m as expected, since the surrogate data is 
supposed to be random in nature, and its D corr should be 
infinite |76| . The saturation of D corr in our case indicates 
that the system dynamics is deterministic in nature. 



In Subsection |3.3| we have reported that the fluctua- 
tions show relaxation type oscillations near the bifurca- 
tion point. 

The time period (T) of these relaxation oscillations in- 
creases with increasing DV near the bifurcation point. 
This eventually results in the vanishing of the limit cycle 
behavior beyond a critical DV (Vh). For larger values 
of DV, the autonomous dynamics exhibit a steady state 
fixed point behavior. Time traces from top to bottom 
in the left panel of Fig [21] depict this period lengthen- 
ing of the oscillatory behavior. A systematic analysis of 
the increment in the period (T), presented in Fig. 21 a) 
[right panel], indicates that the autonomous dynamics 
undergo a exponential slowing down. Consequently, the 
In \ V — Vh\ vs T curve can be fitted by a straight line, 
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FIG. 19: wavelet contour plot of the fluctuations at different different DV shown in Fig |ll| 




FIG. 20: Correlation dimension obtained from multifractal 
spectra (— • — ) and from Grassberger-Procaccia algorithm 
(open circle). 



where Vjj is the bifurcation point separating the oscilla- 
tory domain and the steady state behavior. The results 



of Fig 21 indicate that the system dynamics undergo a 
homoclinic bifurcation [17J at Vh resulting in the loss 
of oscillations. The DV (Vh) divides the plasma into 
two distinct region, The floating potential fluctuation ex- 



hibits relaxation oscillations on the one side of the Vh 
and stable fixed point on the other side, which is termed 
as being in an excitable state and is useful to study noise 
invoked resonance effect in glow discharge plasma that 
will be discussed in Section [4] 



3.4. Conclusions 



Glow discharges are simple systems which exhibit ex- 
otic features depending on the configuration, initial pa- 
rameters, etc. Since various configurations are used in 
different applications like dusty plasma, plasma process- 
ing, etc. it is necessary to understand the plasma dynam- 
ics of these systems as much as possible. However the un- 
derstanding of the complexities in the plasma dynamics 
is quite a challenging job as they arise from many degrees 
of freedom like different sources of free energy, different 
types of wave particle interaction and many other insta- 
bilities. In our present work nonlinear time scries anal- 
ysis has been used to quantify and differentiate complex 
and coherent processes at different parametric conditions. 
We also observed that the chaotic state driven by the DV 
tends to be stabilized through homoclinic bifurcation. 
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4. STOCHASTIC AND COHERENCE 
RESONANCE IN GLOW DISCHARGE PLASMA 

5. INTRODUCTION 

Stochastic resonance (SR) which has been observed in 
many physical, chemical and biological systems [93 [98j 
ESllIDDllIDIllID3IID31IIiaiID3IIDg,isa phenomenon in 
which the response of the nonlinear system to a weak pe- 
riodic input signal is amplified or optimized by the pres- 
ence of a particular level of noise [97], i.e., a previously 
untraceable subthreshold signal applied to a nonlinear 
system, can be detected in the presence of noise. Fur- 
thermore, there exists an optimal level of noise for which 
the most efficient detection takes place. 

In 1981, Benzi and his coworkers, had introduced the 
idea of the SR to model the existence of the 100k year 
cycle of occurrence of the ice age and relatively warm age 
with a temperature difference of the order of 10K, which 
is surprisingly in phase with the Milankovitch cycle. 
They successfully explained the phenomena, assuming a 
bistable system of the global climate, to be separated by 
10K temperature, the Milankovitch cycle as a weak pe- 
riodic forcing and the short term climate fluctuations as 
noise, which has been discussed in Ref E3 EH Ell- 
in 1993, Gang, et al. [TD7] . and in 1994, Kurt Wiesen- 
feld, et al. |108j . had shown that SR is also possible in 
different classes of dynamical systems based not on bista- 
bility but on excitable dynamics. They proposed a sys- 
tem consisting of a potential barrier and above this bar- 
rier the system shows deterministic dynamics (limit cycle 
oscillation) and below, stable fixed point as illustrated in 
Fig [22] 

The figure shows that when the control parameter (po- 
tential) crosses the PB, i.e., the potential is greater than 



PB, the system shows a limit cycle behavior and when 
the potential is less than PB, the system exhibits a stable 
fixed point behavior. Now if we set the control parameter 
(potential) below the barrier and apply noise (stochastic 
perturbation on the potential) to the system and when- 
ever the barrier is crossed, the system returns to its fixed 
point or "rest state" deterministically |107[I108[I109] . i.e., 
whenever the barrier is crossed, the system traverses one 
oscillation. So for moderate noise level we can get back 
relatively coherent oscillations, which is called coherence 
resonance (CR) |110j . Again when a subthreshold pe- 
riodic signal in the form of pulse and noise are added 
to the system below the potential barrier the probabil- 
ity of crossing the barrier by noise at the time of occur- 
rence of the periodic pulse is maximum and hence one 
can get back deterministic dynamics in the form of pe- 
riodic oscillations of frequency of the applied periodic 
pulse for optimum noise level which is the stochastic res- 
onance of an excitable system. So for excitable systems 
both SR and CR are possible. Based upon excitability 
CR have been observed in plasma |111L 1112] , electronic 
circuits [113] , optical systems |114[ I115j . and chemical 
systems [TTSJ EEI3 CH EH] and SR has been observed 
in chemical systems |118j , human brain and many other 
systems [91AIT06]. 

Lin I and Jeng-Mei Liu [111J and A. Dinklage, et 
al. |112j . had studied CR in weakly ionized rf magne- 
toplasma and neon glow discharge plasma based on ex- 
citable dynamics. In both the cases, the excitability has 
been achieved through Hopf bifurcation for particular 
parametric region of the discharge. In our experiments, 
we get excitable dynamics, for certain discharge param- 
eters, in the region greater than the Paschen minimum, 
where the system exhibits selfsustained relaxation oscil- 
lations, and change to excitable fixed point behavior with 
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Potential > V H Limit cycle 



Potential < V H Stable fixed point 



FIG. 22: It shows that when the control parameter crosses the potential barrier (PB), i.e., potential >PB system behaves 
shows limit cycle behavior (up) and when (potential < PB) system shows stable behavior (down). 



increase in the DV, which has been discussed in Subsec- 
tion 13.31 of Section [3j The DV at which these oscillations 
cease may be termed the bifurcation point (Vh) which 
act as PB in this experiment. Near The Vh, the time pe- 
riod of these relaxation oscillations increases dramatically 
upon further incrementing DV. This eventually results in 
the vanishing of the limit cycle behavior beyond the crit- 
ical DV (Vh) and exhibits excitable dynamics through 
homoclinic bifurcation. So below Vh the system shows 
excitable behavior and useful to study noise invoked res- 
onance experiments, in these regions. 



and 5 are the inter-peak interval of the response signal 
and mean peak interval of the subthreshold periodic sig- 
nal respectively. Fig 23 'd) shows that the experimentally 
computed AMD versus noise amplitude D curve has a 
unimodal structure typical for the SR phenomena. The 
minima in this curve corresponds to the optimal noise 
level for which maximum information transfer between 
the input and the output takes place. 



5.1. Stochastic resonance 



5.2. Coherence resonance 



For our experiments on stochastic resonance, the refer- 
ence voltage Vq was chosen such that Vq > Vh and there- 
fore the autonomous dynamics, by virtue of an underly- 
ing homoclinic bifurcation, exhibit steady state behav- 
ior. The discharge voltage V was thereafter perturbed 
V = Vq + S(t) + D£_, where S(t) is the subthreshold pe- 
riodic pulse train chosen for which V = Vq + S(t) > Vh, 
(subthreshold signal does not cause the system to cross 
over to the oscillatory regime) and D£ is the added Gaus- 
sian white noise £ with amplitude D. Subthreshold pe- 
riodic square pulse of width 20 (is and duration 2 ms 
was constructed using Fluke PM5138A function genera- 
tor. Meanwhile, the gaussian noise produced using the 
HP 33120A noise generator was subsequently amplified 
using a noise amplifier. 

Fig. |23|a) — (c) show time series of the system re- 
sponse in the presence of an identical subthreshold sig- 
nal for three different amplitudes of imposed noise. The 
subthreshold periodic pulse train is also plotted, in the 
top most graph of the left panel, for comparison pur- 
poses. Fig. 23 'a) shows that there is little correspon- 
dence between the subthreshold signal and the system 
response for a low noise amplitude. However, there is 
excellent correspondence at an intermediate noise ampli- 
tude [Fig. 23 'b)]. Finally, at higher amplitudes of noise 



the subthreshold signal is lost amidst stochastic fluctua- 
tions of the system response [Fig. 23 [c)). Absolute mean 



For the experiments on coherence resonance DV (Vo) 
was located such that the floating potential fluctuations 
exhibit fixed point behavior. In order to minimize the ef- 
fect of parameter drift, a set-point (Vo) quite far from 
the homoclinic bifurcation (Vh) was chosen. Subse- 
quently, superimposed noise on the discharge voltage was 
increased and the provoked dynamics analyzed. The 
normalized variance (NV) which has been discussed in 
subsection |2.3| of Section ?? was used to quantify the 
extent of induced regularity. It is defined as NV — 
std(t p ) /mean(t p ) , where t p is the time elapsed between 
successive peaks. It is evident that more regular the in- 
duced dynamics the lower the value of the computed NV. 
For purely periodic dynamics the NV goes to zero. 



Fig. 24 a) — (c) (left panel) show the time series of 



the floating potential fluctuations for different noise lev- 
els and Fig |24|d) (right panel) is the experimental NV 
curve as a function of noise amplitude D. The point (a) in 
Fig 24 'd) (time series shown in Fig. |24[a)) is associated 
with a low level of noise where the activation threshold is 
seldom crossed, generating a sparsely populated irregular 
spike sequence. As the noise amplitude is increased, the 



NV decreases, reaching a minimum (b) in Fig 24 'd) (time 



series shown in Fig. 24 'b)) corresponding to an optimum 



difference (AMD) which has been discussed in Subsec- 
tion |2.3| of Section [2j used to quantify the information 
transfer between the subthreshold signal and the system 

response, is defined as AMD — abs(mean( t -j — 1)). t p a consequence of the dynamics being dominated by noise. 



noise level where maximum regularity of the generated 
spike sequence is observed. As the amplitude of superim- 
posed noise is increased further, the observed regularity 
is destroyed manifested by an increase in the NV; label 
(c) in Fig 24 'd) (time series shown in Fig. 24 'c)). This is 
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FIG. 23: Noise induced system response in the floating potential fluctuations for low, medium, and high amplitude noise 
in conjunction with a subthreshold periodic square pulse. The right panel [Fig 23 shows the AMD as a function of noise 
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FIG. 24: Emergence of coherence resonance: The right panel shows the NV as a function of noise amplitude for the experiments 
performed at Vo = 344 V and pressure = 0.5 mbar. Left panel: The time series of the floating potential fluctuations for (a) low 
noise, (b) optimum noise and (c) high level noise 



5.3. Discussion 

The effect of noise has been studied experimentally 
near a homoclinic bifurcation in glow discharge plasma 
system |17j . Our study demonstrates the emergence of 
SR for periodic subthreshold square pulse signals and 
the induction of CR via purely stochastic fluctuations. 
In SR experiments, the efficiency of information trans- 
fer was quantified using AMD instead of the power norm 
which has been utilized elsewhere |101j . The advan- 
tage of using this method in comparison to the power 



norm (Co(0)) [101 1 lies in the fact that AMD remains 
independent of the lag between the measured floating 
potential and the applied periodic square pulse. This is 
of relevance to our experimental system, where invari- 
ably there exists a lag, at times varying in time due to 
the parameter drifts. Comparison between the estimated 



C*o(0) and AMD have been shown in Fig 26 It is obvious 
from Fig 26 [ai) that Cn(0) does not show any peak at 
optimum noise level. Whereas, AMD shows nice agree- 
ment with regular spiking of the signal [Fig 26'i>i)]. For 
the CR experiments it was occasionally observed that 
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FIG. 25: Emergence of coherence resonance: NV as a function 
of noise amplitude for the experiments performed, at Vo = 
344 V and pressure = 0.5 mbar. 



while with an initial increase in noise amplitude (D) NV 
reaches a minimum, the subsequent rise of NV for even 
higher amplitudes of noise was suppressed as shown in 
Fig [25] This leads to the modification of the unimodal 
profile, a signature of the CR phenomenon. 

A possible explanation for this suppression is that by 
virtue of the superimposed high frequency noise (band- 
width 500 KHz) and fast responding internal plasma dy- 
namics, the system has the capability of exciting high 
frequency regular modes within the ion plasma frequency 
(105 kHz). This in turn leads to the persistence of low 
NV values. Finally, in Refs [ITTJ [TT2] both the destruc- 
tive and constructive role of noise (CR only) have been 
reported for glow discharge and magnetized rf discharge 
plasma systems respectively. However, both these ex- 
periments were carried out in the vicinity of the Hopf 
bifurcation. In contrast, for the present work we studied 
both SR and CR in the neighborhood of the homoclinic 
bifurcation. 



6. REALIZATION OF SOC BEHAVIOR IN 
GLOW DISCHARGE PLASMA 

7. INTRODUCTION 

Many empirical analyses suggest that power-law be- 
havior in the distribution of some quantities are quite 
frequent in nature for example, 1// noise, Gutenberg- 
Richter law for the magnitudes of earth quakes, stock 
market, forest fires, ecology, weather etc [1201 11211 11221 
11231 1124] . From the point of equilibrium statistical me- 
chanics, power law is considered as a signature of critical 
fluctuations observed near a second order phase transi- 
tion. As we approach the critical temperature, correla- 



tion length diverges and measurable quantities exhibit 
power law. However, power laws are often seen in many 
other systems, not necessarily connected with equilib- 
rium situations, some of these systems are driven sys- 
tems which are very different from equilibrium systems. 
In the case of second order transition, the critical state 
is reached by only through parameter tuning. The idea 
of self organized criticality (SOC) introduced by Bak, 
Tang and Wiesenfeld |125j was that when the systems 
are driven away from equilibrium, they reach a state near 
to, but not at, the state that is marginal to major disrup- 
tions without any significant tuning of the system from 
outside. These systems are characterized by a spectrum 
of spatial and temporal scales of the disruption events 
that exists in remarkably similar form in a wide variety 
of different physical systems. 

The simple concept of SOC is the sandpile. Consider, 
the process of building a sandpile by dropping sand on 
a particular spot. As we add sand, the pile will grow 
and always chooses its shape which is independent of the 
builder. This property is called self organization. The 
sandpile surface will steepen until its slope reaches the 
angle of repose, the critical angle beyond which further 
addition of sand will cause cascades of sand to topple 
down the sandpile. Most parts of the sand pile have a 
slope very close to this critical angle. Any further addi- 
tion of the sand grains will trigger cascades of sand parti- 
cles which are more commonly known as avalanches. The 
sand pile is a good example of a critical state, where the 
correlation length of perturbations extends over the en- 
tire sandpile, i.e., independent of the size of the pile, an 
addition of one grain of sand close to the top of the pile 
can influence grains at the bottom. In this process there 
is no external parameter needs to be tuned because the 
critical state in the sand pile is attained as a consequence 
of the slow addition of the sand grains and their spatial 
redistribution by avalanches |125l I126| . 

The SOC concept has been quite rigorously deployed 
to explain some of the turbulent transport observations in 
magnetically confined fusion devices like Tokamaks. If we 
consider free energy source in the plasma is the gradient 
of the density. When this gradient becomes greater than 
some critical value depending upon plasma parameters, 
by fueling the plasma and turbulence turbulence develops 
around the location. This eddy flatten the density of 
that position and free energy moves to its adjacent place 
and in the similar way eddies may be triggered the new 
location. In this way the density profile is modified again 
and the process continues and produces avalanche-like 
phenomenon which has been successfully explained by 
SOC concept. 

Though SOC concept has been applied quite success- 
fully in high temperature plasmas like tokamak |127[ 1128} 
I129j , we first observed that the low temperature glow dis- 
charge plasmas also exhibits SOC behavior, for the region 
less than the Paschen minimum which has been described 
in Subsection 12.11 of Section 

From the the Paschen curve [Fig [4] of Section|] for our 
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system we have seen that, at low pressure i.e., for less 
the than Paschen minimum, discharge struck at higher 
discharge voltage (DV) [subsection 2.1 of Section In 
the lower region, keeping the DV at a constant value of 
800 V, a discharge was struck with a very faint glow at 
~ 9 x 10~ 3 Torr, and then as the filling pressure was grad- 
ually increased by means of the needle valve, the intensity 
of the glow suddenly became bright at about 1.6 x 10~ 2 
Torr. The sudden change in the plasma dynamics at the 
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FIG. 27: The variation of plasma discharge current (solid 
line), and plasma floating potential (dotted line) with pres- 
sure. 

critical pressure led to a change in the floating potential, 
and the discharge current. The variation of the plasma 
floating potential, and the plasma discharge current with 



pressure have been shown in Figure 27 It shows that, the 
current (solid line) jumped to a larger value by a factor 
of 15 at 1.6 x 10 -2 Torr, and then increased gradually 
with pressure. On the other hand, the floating poten- 
tial rapidly fell to a negative value (« —9 V) up to the 
same critical pressure (1.6 x 10~ 2 Torr), and then again 
increased (dotted line in Fig. 27) with further increase in 
pressure, until it finally settled down to « —5 V. So the 
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FIG. 28: Plot of critical pressures with discharge voltages. 
Corresponding floating potential, and discharge are within 
brackets. The first value in the brackets is the floating poten- 
tial, and the second one is the discharge current in (V, mA) 
unit. 



critical pressure divides the scanned region into two re- 
gions I and II, shown by a vertical line (— .— ) in Figure 27 
The plasma density, temperature, and electron-electron 
collision mean free path in region II are 10 7 — 10 8 cm~ 3 , 
2—4 eV, and 2.6 x 10 5 cm respectively. However, in region 
I, it was almost impossible to obtain the I-V characteris- 
tics because of their extremely low values. Qualitatively, 
looking at discharge current, and glow intensity, region II 
is probably a normal glow discharge region, while region I 
might be the dark or subnormal glow discharge region. In 
the region I system exhibits SOC behavior. The bound- 
ary between the two region (critical pressure) where the 
transition takes place is not a fixed point, but decreases 

It is 



with the discharge voltage as seen in Figure 28 



also seen that the floating potential, and the discharge 
current also decrease with discharge voltage. 

Figures 29 a) , (b) , and (c) are the typical electro- 
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FIG. 29: Figure shows the electrostatic fluctuations at pres- 
sures 0.9 x 1(T 2 (a), 1.2 x 1CT 2 (b), 1.5 x 1CT 2 (c), 1.7 x 1CT 2 
(d), and 2.2 x 10" 2 Torr (e) respectively. 




FIG. 30: ln(power) vs. In / plot. The solid line shows the 
best fit. 



static fluctuations at pressures 0.9 x 1CP 2 , 1.2 x 1CP 2 , 
and 1.5 x 1CP 2 Torr respectively in region I, while (d), 
and (e) are the fluctuations at pressures 1.7 x 1CT 2 and 
2.2 x 1CP 2 Torr respectively in region II. In the next sec- 
tion it will be shown using different analysis techniques 
that the fluctuation in the region I is consistence with 
SOC behavior. 



7.1. Analysis of SOC behavior 

The experimental evidences considered as main ingre- 
dients of SOC are l//' 5 (/3 >0) power law (where / 
is the frequency of the fluctuations obtained from fast 
Fourier transform) |125l 11301 11311 1132) , long-range cor- 
relation [133] . and nongaussian probability distribution 
function (PDF) [134] . From the power spectral analysis 
which have been discussed in Subsection |2. 3| of Section [2] 
we have estimated the (3 from ln(Power) versus ln(/). 

For long-range time correlation we estimated the Hurst 
exponent H using rescale range method, and the expo- 
nent (a) of autocorrelation function (ACF) decay, as de- 
scribed below. Both the techniques have been discussed 
in Subsection 12.31 of Section |2] 

In order to verify nongaussianity we obtained the PDF 
of the fluctuating data. 



Figure 30 shows In (Power) vs. In / of the fluctuations 
in region I, from which we estimated the exponent to be 
w 1.6. This agrees well with the numerical |125|,I135] . as 
well as experimental observations |1301 11311 1133j in the 
presence of SOC behavior. 

31 shows a typical plot of ln(R/S) vs. ln(timc 
ic fluctuations in region I, for a pressure of 
The Hurst exponent H is about 0.96 ± 0.01 



Figure 
lag) of t 
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FIG. 31: ln(R/S) vs. ln(lag) plot for the electrostatic fluctu- 
ation at the pressure 1.4 x 10~ 2 Torr. The solid line shows 
best fit. 



gion I. On the other hand, for sinusoidal like oscillation 
in region II, the Hurst exponent is 1 for the lag length 
of one period of oscillations, and zero with more than 
one period lag |136j as shown in Figure 32 
same region (II) multi slope ln(R/S) vs. 



Also in the 
ln(lag) plot as 
seen in Figure [33] probably due to plasma instabilities 
of many frequencies |136j . The ACF exponent a has 
been calculated for the fluctuations in region I, from the 
In (ACF) vs. In (time lag) plot as shown in Figure I 



34 The 



(this indicates long-range time correlations) |127j , and it 
is almost constant over the entire pressure range in re- 



ACF (Fig. |34| ) shows power law up to about 6 times the 
decorrclation time, and after that it follows exponential 
decay. Average value of a is about 0.30. Using the re- 
lation H = (2- a)/2 [I57j . H calculated from ACF is 
~ 0.85, which is close to the value of H, calculated using 
R/S technique. In region II, no power law decay of ACF 
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In(lag) 

FIG. 32: R/S as a function of time lag for the electrostatic 
fluctuation of the coherent oscillations at the pressure 1.7 x 
1CP 2 Torr. For one time period lag, H=l, and for lag more 
than one period, H=0. 




In(lag) 

FIG. 34: ln-ln plot of ACF vs. time lag at pressure 1.4 x 1CT 2 
Torr. Up to 6 decorrelation times, it shows power law (dotted 
line), and after that it follows exponential decay. 




been shown in Table III for pressures 0.9 x 10 2 , 
10~ 2 . 1.5 x 10~ 2 Torr. 



1.2 x 



In(lag) 

FIG. 33: ln(R/S) versus In (lag) plot at the pressure 2.2 x 10~ 2 
Torr, more than one slope indicates instabilities with many 
modes. 



has been observed. The PDF of the floating potential 
fluctuations in region I, seen in Fig. 35 'a) clearly shows 
a nongaussian nature. Corresponding best gaussian fit is 
given by dotted curve in the same figure. 

We suspect that there might be a slightly bimodal dis- 
tribution similar to Ref . |131] . Figure 35 'b) shows the 



gaussian nature of the fluctuations in region II. 

Our results of Hurst exponent, H > 0.5, ACF expo- 
nent, a ~ 0.30, nongaussian PDF, and power spectral 
index f3 ~ 1.60 in the pressure range 9 x 10~ 3 — 1.6 x 10~ 2 
Torr, are consistent with the systems exhibiting SOC like 
behavior. 

Comparison of a, H by ACF, H by R/S, and (3 have 



7.2. Conclusion 

We have obtained SOC behavior over a finite range of 
neutral pressure of 9 x 10~ 3 — 1.6 x 10~ 2 Torr for a fixed 
discharge voltage and this region belongs to the region 
less than the Paschen minimum |19j . Glow discharges 
are simple systems, but their physics can be quite com- 
plicated due to the presence of several phenomena like 
avalanche breakdown, ionization waves, low frequency 
ion-acoustic instability, double layer, chaos, etc. Most 
of them are highly nonlinear processes and hence one re- 
quires different techniques both statistical and spectral 
to investigate and understand their behavior. From our 
present analysis we observe that the plasma dynamics in 
the region I, is compatible with self organized criticality, 
while region II, is not. This could also imply that plasma 
transport in region I, is quite different from region II. 



TABLE III: In the following table the decay exponent a of 
the ACF, H from ACF, H using R/S, and the power spectral 
index j3 have been shown in the second, third, fourth, and 
fifth column respectively, for the pressures shown in the first 
column. 

pressure a Hurst Hurst f3 
10- 2 (Torr) (ACF) (ACF) (R/S) (PS) 



0.9 


0.38 


0.81 


0.97 


1.85 


1.2 


0.34 


0.83 


0.96 


1.65 


1.5 


0.23 


0.88 


0.95 


1.60 
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FIG. 35: The probability distribution function (PDF) of the 
fluctuation at pressure 1.2 x 10 -2 Torr, and the dotted line is 
the corresponding gaussian fit. 



8. SUMMARY OF THE WORK 

8.1. Summary 

The thesis presents the experimental observations of 
the nonlinear phenomena— chaos to order transition, self 
organized criticality behavior and noise invoked processes 
like the stochastic and coherence resonance in an argon 
dc glow discharge plasma [Fig l]of Section[2]. Depending 
upon the gas pressure (p), a discharge was struck at dif- 
ferent discharge voltages (DV) , which exhibits a Paschen 
curve like behavior [Fig. [4] of Section [2] . The region left 
to the Paschen minimum is narrower than the region to 
the right and the behavior of the plasma floating poten- 
tial fluctuations were different on these two sides of the 
Paschen curve. 

When the system was operated for the pressures 
greater than the Paschen minimum, an anode glow or an- 
ode spot was observed to form on the anode and with the 
increase in the DV its size decreased and finally dimin- 
ished at a certain DV. The annular radii of the glow esti- 
mated using CCD camera was in the range of 1.3 — 0.32 
mm. This observation is very consistent with the the- 
oretical estimations of thickness (d) using the relation 



3.7 x 10 



-6 kT 

iP' 



where k, T, <Ji, and P are the Boltz- 



mann constant, room temperature in Kelvin scale, ioniza- 
tion cross-section and pressure in mbar respectively. The 
estimated thickness of the anode glow using the above 
relation for P = 0.95 mbar, T — 300 K, is approximately 
0.81 mm is within the range of the thickness as obtained 
from the experiment (1.3 — 0.32 mm). 

An interesting feature associated with the anode glow 
was the different types of oscillations in the floating po- 
tential at different pressures, and these fluctuations have 
been analyzed for the three typical pressures in the re- 
gion greater than the Paschen minimum. At about 0.89 



mbar (pd « 20.02 mbar-mm), the discharge was initi- 
ated at w 288 V and simultaneously irregular relaxation 
type of oscillations were observed and the regularity of 
the oscillation also increased with the DVs. Around 509 
V, these oscillations disappeared and this point (DV) is 
termed as the bifurcation point (Vh)- Almost similar 
behaviour had been observed at 0.95 mbar (pd = 21.37 
mbar mm) and at 1 mbar (pd =22.5 mbar mm ). Ini- 
tially, the power spectrum of the oscillations as of broad 
band nature indicating chaotic behavior of the system. 
The frequency of the oscillations were around the ion 
plasma frequency. An estimate of the frequency of these 
instabilities can be obtained fr om th e ion transit time in 
the plasma r(d) = y^— = dj y^p^ where d is the elec- 
trode distance. The estimated ion transit frequency (-) 
for our experimental system is ~ 19 kHz which agrees 
well with the frequency of the relaxation oscillations of 
the floating potential. 

The presence of relaxation oscillations has been at- 
tributed to the formation of highly nonlinear structures 
like double layers. We therefore tried to estimate the 
correlation dimension (D corr ) and the +ve Lyapunov ex- 
ponent (Al) of all the signals. The estimated D corr for 
the data for the three pressures (0.89, 0.95 and 1.0 mbar) 
were greater than 3.8 to begin with and decreased with 
increase in DV. Since D corr is a measure of the complex- 
ity of the system, initially, the system was in a complex 
state as D corr for all the three pressures were high and 
decreased in complexity with increase in DV. 

The presence of a +ve Lyapunov exponent (Xl) is the 
most reliable signature of the chaotic dynamics. The 
positive Al has been identified for DV 283, 284, and 290 
V at 0.95 mbar and for 293, 296, 300 and 305 V at 1 mbar 
respectively. At higher pressures we find that in general 
Al is +ve and D corr > 3, suggesting a low dimensional 
chaos. Wavelet analysis also showed presence of chaos at 
the initial stage of the discharge. 

Generally, the floating potential fluctuations were com- 
plex in nature at the initial stage of discharge and became 
regular with increase in DV and converted to relaxation 
oscillations and vanished through homoclinic bifurcation 
at Vh- Applying noise beyond the Vh noise invoked ex- 
periments, coherence and stochastic resonance had been 
performed. For the coherence resonance we had only ap- 
plied noise to the system to get back the autonomous 
dynamics, whereas, for the stochastic resonance, noise 
and periodic pulse had been applied simultaneously. In 
order to overcome the plasma related problem in analyz- 
ing the stochastic resonance, we had devised new analysis 
techniques absolute mean difference. 

When the system was operated in the pressure region 
less than the Paschen minimum, for small range of p 
(0.9 - 1.5 x 10~ 2 Torr or 1.2 - 2 x 10~ 2 mbar), it was 
observed that the behavior of the floating potential fluc- 
tuations was consistent with SOC. In order to establish 
the SOC behavior, we had checked the power law be- 
havior of the power spectrum, the presence of the long- 
range correlation by estimating Hurst exponent (H) us- 
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ing R/S technique and the exponent (a) of ACF decay, 
and the nongaussian probability distribution function. 
The results of Hurst exponent (0.96 ± 0.01) greater than 
0.5, ACF exponent, a ~ 0.30, nongaussian PDF, and 
power spectral index /?wl.7±0.1in the pressure range 
9x 10~ 3 — 1.6x 10~ 2 Torr, are consistent with the systems 
exhibiting SOC like behavior. 

8.2. Scope of the future works 

Though we have investigated some of the nonlinear 
processes during the the course of this thesis work, there 
are, many problems which we plan to take up as future 
plans. They are as follows. 

• Effect of the noise on the autonomous dynamics. 

• Effect of suprathreshold signal and noise to the plasma. 

• Investigation of the existence of canard orbits in 
plasma. 

• Non-chaotic attractors may be investigated by apply- 
ing two non-commensurate periodic signals. 

• Chaos control and synchronization. 

• Modeling of the experimental results. 

Dusty plasma which can be produced very easily in 
our experimental system, is another area where, a lot of 
nonlinear phenomena can be explored. 
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